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Abstract 

All  systems  are  subject  to  failure  and  must  be  repaired  to  be  kept  in  work¬ 
ing  order.  The  word  'repair'  is  used  here  in  a  broad  sense.  It  can  consist  of 
replacement  with  a  brand  new  item,  or  checkups  at  periodic  intervals,  or  several 
other  forms  of  'minimal'  repairs.  In  this  paper,  we  describe  several  kinds  of 
'r-pilrs'  and  maintenance  models.  After  studying  some  properties  of  the  sto¬ 
chastic  process  of  failure  times,  we  compare  different  maintenance  models  by 
comparing  the  expected  number  of  failures  in  time  t. 
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1.  Introduction. 


All  systems  are  subject  to  failure  and  must  be  repaired  to  be  kept  in  work¬ 
ing  order.  The  word  'repair*  is  used  here  in  a  broad  sense.  It  can  consist  of 
replacement  with  a  brand  new  item,  or  checkups  at  periodic  intervals,  or  several 
other  forms  of  'minimal'  repairs.  In  this  paper,  we  describe  several  kinds  of 
'repairs'  and  study  some  properties  of  the  stochastic  process  of  failure  times 
under  various  maintenance  models. 

Consider  a  system  consisting  of  only  one  unit.  This  unit  is  put  into  use 
at  time  t  =  0  and  has  life  distribution  P.  IVhen  it  fails,  we  can  perform  a 
'perfect'  I'spair  or  a  'minimal'  repair.  'Perfect'  repair  means  that  a  new  unit 
’  '  -  '  ■  'o  dis r ’bui-'' on  is  ^  is  put  in  the  place  of  the  failed  unit.  A  'G-mini- 

Tir  ’  “  pair'  means  that  the  failed  unit  is  replaced  by  a  unit  with  life  distribu¬ 
tion  "  and  ago  equr.Js  to  t!.e  '''ffsetive  age'  of  the  failed  unit.  More  formally, 
let  the  failure  time  of  the  x'ailed  unit  as  computed  from  time  zero  or  the  last 
perfect  repair,  whichever  came  last,  be  y.  Then  the  probability  that  the  life 
of  the  G-minimal  repaired  replacement  exceeds  x  is  G(x  +  y)/G(y)  for  xSO.  This 
definition  of  minimal  repair  was  proposed  by  Ascher  (1978)  and  has  been  used  by 
Brown  and  Proschan  (1983). 

Brc'.vn  and  Proschan  (1983)  considered  the  following  maintenance  model.  A 
-rerri-ting  of  a  single  unit  starts  out  with  a  unit  whose  life  distribu- 
t'.'  I  .  V.hcnever  a  failure  occurs,  a  coin  with  probability  p  for  heads  is 

tv.  i..d2fcr.dcnt  of  previous  history.  If  the  coin  comes  up  heads  a  perfect 
repair  ^ ’'rrformed.  Otherwise  an  F-minimal  repair  is  performed.  Notice  that 
each  cp.-.ch  perfect  repair  is  a  rcpc/’.c.v  tion  point  for  the  process  of  failure 
times.  lr.'vn  and  Proschan  (1983)  obtained  the  distribution  of  the  time  between 


perfect  repairs  (which  is  also  the  same  as  the  waiting  time  for  the  first  per¬ 
fect  repair)  and  established  some  of  its  monotonicity  properties  in  terms  of 
similar  properties  of  F. 

It  is  easy  to  see  that  the  distribution  of  the  time  between  perfect  repairs 
depends  in  a  simple  fashion  on  the  distribution  of  the  ninnber  of  minimal  repairs 
before  time  t  in  another  maintenance  model  called  the  Forever  Minimal  Repair 
(Ff-IR)  model  where  only  G-minimal  repairs  are  performed  (see  Theorem  2.1).  This 
leads  to  an  alternative  derivation  of  Lemma  2.1  of  Brown  and  Proschan  (1983). 

What  are  the  interesting  questions  concerning  maintenance  models?  By  very 
definition,  a  system  can  be  maintained  indefinitely  under  all  repair  models,  if 
the  lifetimes  of  repaired  units  are  unbounded  to  the  right.  Some  models  will 
be  more  expensive  to  maintain  than  others.  In  this  paper  we  will  study  the 
failure  time  processes  of  some  maintenance  models  and  compare  the  expected  ntJm- 
ber  of  failures  in  tree-  t,  which  is  roughly  related  to  the  cost  of  maintaining 
the  system  till  time  t. 

A  maintenance  model  where  perfect  repairs  are  performed  at  each  failure  is 
called  the  Her.ove  and  Replace  (RAR)  model.  The  failure  times  in  this  model  are 
easily  studied  by  using  standard  renewal  theory.  In  Section  2  we  compare  the 
expected  number  of  failures  in  the  RAR  and  FMR  models.  It  is  shown  that  the  RAR 
model  has  smaller  expected  number  of  failures  if  the  unit  life  distribution  is 
IrRA.  '3n  the  other  hand,  the  FMR  model  has  a  smaller  expected  number  of  failures 
if  the  unit  life  distribution  is  DFRA. 

In  Section  3  we  consider  other  maintenance  models  where  two  types  of  mini¬ 
mal  repairs  are  available.  These  minimal  repairs  may  be  chosen  at  random  before 
the  start  of  the  process  and  fixed  ever  after  (one-shot  random  repair  model),  or 
can  be  chosen  at  random  at  start  and  alternated  thereafter  (alternating  repair 


model)  or  can  be  chosen  at  random  at  start  and  again,  independently,  at  each 
failure  (completely  random  repair  model) .  We  then  compare  the  expected  number 
of  failures  in  time  t  under  all  these  maintenance  models. 


2.  Maintenance  models  with  only  one  type  of  minimal  repair. 

In  this  section  we  consider  maintenance  models  with  only  one  type  of  mini¬ 
mal  repair.  Of  course,  if  we  do  not  wish  to  use  minimal  repair,  we  can  use  per- 
feet  repair,  which  is  assumed  to  be  always  available. 

Consider  the  Remove  and  Replace  (RAR)  model  in  which  a  failed  unit  is 
replaf o(!  with  a  new  unit;  that  is  a  perfect  repair  is  performed  at  each  failure 
The  lifetimes  ox  all  the  units  are  independent  and  identically  distributed  with 
common  distribution  function  F.  Then  the  ntimber  of  failures  before  t,  N(t),  is 
a  standard  yeriowal  process.  It  is  vjell  known  that  the  probability  generating 
function  t|i(t,  q)  of  K(t)  is  given  by 

00 

Kt.  q)  =  Fj{t)  +  I  q"[F^^i(t)  -  F^(t)] 
n  =  1 

09 

=  (l-q)  I  q"‘X(t), 
n=  1  " 

— *  #00 

where  F  denotes  the  n-fold  convolution  of  F.  Furthermore,  if  J  F(t)dt = p < “, 
n 

then  ti;e  renewal  theorem  states  that 


We  refer  to  Barlow  and  Proschan  (1975),  p.  167  for  the  proof  of  these  results. 
These  well  known  results  will  be  used  as  yardsticks  for  later  comparisons. 

In  contrast  to  the  above  model  consider  the  Forever  Minimal  Repair  (FMR) 
model  in  which  each  failed  unit  is  G-mininally  repaired.  More  formally,  if 
there  is  a  failure  at  time  t,  then  it  is  replaced  with  a  working  unit  with  life 
distribution  G  and  age  t,  i.e.,  with  a  random  life  T  whose  distribution  is  given 
by 

P(T  >  s)  =51^-111  ,  S2  0. 

G(t) 

W3  also  ass’.ime  that  the  life  distribution  of  the  first  unit  is  G.  Let  N(t)  be 
the  number  of  failures  before  t.  Theorem  2.1  below  obtains  the  probability 
generating  function  cf  N(t). 

Theorem  2.1.  The  probability  generating  function  of  N(t)  is  given  by 
(2.3)  <>(t,  q)=GP(t)  where  p=l-q. 


Proof. 


th 


Let  X  be  the  time  of  the  n  failure  and  let  F  be  the  distribution  of 
n  n 

X  .  Then  for  r  >  1 , 
n 


(2.4) 


P[N(t)  =n]  =F^^jCt)  -F^(t) 


=  /;nx„.i2t|X^  =  u)dF„(u) 


Define  a^^Ct)  =  1  and 


t 

-  for  n  Si. 

"  "  G(u) 


Then  P[N(t)  =  n]  =  a^(t)G(t)  for  nsO.  Differentiating  (2.4)  and  (2 
respect  to  t,  we  obtain 


-dF 


n+1 


(t)  +dF^(t)  =G(t) 


G(t) 


-  a^(t)dG(t). 


dF  ,  (u) 
G(u) 

a^(u)dG(u) 

G(u) 


Let  A(t,  q)  =  \  q”ci  (t) .  Then 

n  =  0 


A(t,  q)  =  a^(t)  + 


0  I  ''"ViW 

n  =  0 


=  l  +  q/^ACu,  q) 

^  G(u) 


Differentiating  \Jith  respect  to  t,  this  yields 

A'(t,  q)=qA(t,  q)^-— -  • 
G(t) 

Since  A(0,  q)  =  1,  this  has  a  unique  solution  given  by 

A(t.  q)  =G’‘’  (t). 


5)  with 


i 


-  6  - 


Hence  q)  =  I  q’'p[N(t)=nl 
n  =  0 


00 


=  I  q"G(t)a  (t) 
n  =  0  " 


=  G(t)A(t,  q) 

=  GP(t).  n 

Thus  in  the  FMR  model,  the  expected  number  of  failures  before  t  is  given  by 

(2.6)  E(N(t))  q)lq^  i  =  -JnG(t). 

Using  (2.2)  and  (2.6),  we  can  compare  the  expected  number  of  failures  before 
t  in  the  RAR  model  and  FMR  model.  The  expected  number  of  failures  before  t  is 
smaller  in  the  RAR  model  (for  large  t)  if  and  only  if 

(2.7)  G(t)  ^ 

Let  G  be  an  IFRA  distribution  with  mean  u.  Then  -^logG(t)  is  increasing  and 
cannot  lie  always  below  l/p  on  always  above  l/;i,  because  this  would  contradict  the 
assumption  that  the  mean  of  G  is  u.  Thus  (2.7)  is  true  for  all  large  t  and  the 
expected  number  of  failures  before  t  is  smaller  in  the  RAR  model  than  in  the 
Ff!R  model  if  G  is  IFRA.  Similarly,  if  G  is  DFRA  then  the  expected  number  of 
failures  before  t  is  smaller  in  the  B'R  model  than  in  the  RAR  model,  for  all 
large  t, 

VJe  can  add  a  slight  element  of  generality  to  the  Ff-1R  model  by  assuming  that 
the  first  unit  has  life  distribution  F  and  that  G-minimal  repairs  are  performed 
at  all  failures.  I'^e  will  call  this  the  Extended  Forever  Minimal  Repair  (EFMR) 
model.  The  probability  generating  function  of  N*(t),  the  number  of  failures 
before  t  in  this  model  is  obtained  below. 


Corollary  2.2.  The  probability  generating  function  of  N*(t)  is  given  by 


«*(t. 


Proof. 

OD 

Since  (J>*(t,  q)  =  I  PtN*(t)=n]q" 
n  =  0 

00 

=  FCt)^  I  q"/!  P[N(t- s)  =n-  l]dF(s) 
n=  1  ^ 

=  F(t)  +q/J^'(t  -  s,  q)dF(s). 

where  <i(u,  q)  and  N(u)  are  as  given  in  Theorem  2.1  except  that  the  first  unit  has 
life  distribution  G  with  age  s,  and  when  a  failure  occurs  at  time  w,  the  failed 
unit  is  replaced  by  a  unit  with  life  distribution  G  and  age  w  +  s.  From  (2.3), 
wc  have 

^*(t,q)  =F(t)  dF(s).  P 

°^(S) 

If  F  and  G  have  proportional  hazard  rates,  i.e.  G{t)  =F^(t)  for  some  y>0, 
then  the  above  probability  generating  function  becomes 

«*(t,  q)  =  (l-Y^)F(t)^j-3^F^P(t). 

Though  this  is  a  distribution  function  in  t  for  all  y  >  <^ ,  it  is  a  convex  combi¬ 
nation  of  F  ard  f'*^^  only  for  0<y<1. 

The  stochastic  processes  (NCt),  t  ^  0}  and  {N*(t)>0}  are  interesting  pro¬ 
cesses  in  their  own  right,  n-iorrer.  2.1  shows  that  the  marginal  distribution  of 
N(t)  is  Poisson  with  parameter  -inG(t).  The  number  of  failures  before  time  t, 
N(t),  in  the  FMR  model  is  exactly  the  number  of  record  values  before  time  t  in  a 


sequence  of  i.i.d,  random  variables  with  common  distribution  G.  More  general 
results  about  Nft)  can  be  found  in  the  literature  for  record  values.  For  example, 
ShoTTOck  (1972)  has  shown  that  {N(t),  t > 0}  is  a  non-homogeneous  Poisson  process. 
Our  proof  of  Theorem  2.1  can  be  modified  to  establish  the  same  result. 

We  can  interpret  the  probability  generating  functions  obtained  in  the  P1R 
and  EFMR  models  as  the  distribution  function  of  the  waiting  time  between  two 
perfect  repairs  in  another  maintenance  model  due  to  Brown  and  Proschan.  In  the 
Brown  and  Proschan  model,  the  system  starts  out  with  a  unit  with  life  distribu¬ 
tion  F.  At  each  failure,  a  coin  with  probability  p  for  heads  is  tossed  independ¬ 
ent  of  previous  history.  If  it  is  a  head,  a  perfect  repair  is  performed  and  the 
failed  unit  is  replaced  by  a  unit  whose  life  distribution  is  F.  If  the  coin 
turns  up  tails,  tlie  unit  is  replaced  by  a  unit  with  life  distribution  G  and  age 
t*  where  t*  is  the  time  from  the  beginning  or  the  previous  perfect  repair,  which¬ 
ever  came  last.  In  short,  when  the  coin  turns  up  tails,  a  G-minimal  repair  is 
performed.  !‘'c  distinguish  two  models,  the  Brown -Pros chan  (BP)  model  in  which 
F  =G  and  the  Extended  Brown-Prcschan  (EBP)  model  in  which  F  and  G  are  not  equal. 

In  both  models,  the  epochs  of  perfect  repair  form  regeneration  points  and  the 
process  starts  over  again.  Thus  it  vjould  be  interesting  to  obtain  the  distribu¬ 
tion  of  the  waiting  time  for  the  first  perfect  repair  or  the  time  between  two 
successive  perfect  repairs.  Denote  this  waiting  time  by  U  in  the  BP  model  and 
bv  U'*  in  the  EBP  r;odel. 

(i)  P(U>t)  =<>(t,  q)  =GP(t). 


Theoren  2.3. 


9  - 


Proof. 

prove  only  (i) .  The  proof  of  (ii)  is  similar. 

00 

Since  P(U>t)  =  5!  P(exactly  n  failures  occurred  before 

n  =  0  t  and  a  G-minimal  repair  is  performed 
at  each  failure) 

00 

=  I  q"P[N(t)=n], 
n  =  0 


where  N(t)  is  the  number  of  failures  before  t  in  the  FMH  model.  From  the  defi¬ 
nition  of  <{i(t,  q) ,  we  have 


p(u>t)  =  <>(t,  q).  n 

This  provides  an  alternative  proof  for  Lemma  2.1  of  Brown  and  Proschan 
C1983)  . 

Brown  and  Proschan  studied  the  ageing  property  of  the  distribution  of  U  in 
terms  of  those  of  G.  \)e  give  a  result  below  that  gives  an  upperbound  for  EU* 
a.'^ri'ming  that  G  is  N3U  or  OMRL. 


Theorem  2.4.  Let  G  be  N’h'J  or  DMRL.  Then  under  the  EBP  model, 

EU*  <  u  +  for  0  <  p  <  1 
P  ' 

where  u  =  jF(tlet  and  ri  =  jG(t)dt. 

Fr "of . 

It  follovfs  from  Brown  and  Proschan  (1983)  that  is  NBUE  and  /G^(t)dt  ^ 


(1/p)  /  G(t)dt  for  fl  <  p  <  1 .  Thus 


s 


5  in/p  if  u  2  0  and  0  <  p  5  1 


By  Theorem  2.3,  we  have 


EU*  =  /4.*(t,  q)dt 

°  ‘^gP(u) 

=  V  ♦  q /^  [ /"  ^  (O  dt/gP  (u)  ]  dF  (u) 
Sw  +  q/QGP(t)dt 


5  V  +-^111 .  0 

P 


3.  Mainterance  models  with  two  tvnes  of  minimal  renairs. 


We  extend  the  results  of  the  previous  section  to  models  where  two  types  of 
minimal  repairs  are  available.  This  could  correspond  to  the  practical  situation 
where  repairmen  with  two  types  of  training  are  called  in  to  do  the  repair,  or 
where  repairs  are  performed  with  spares  from  deteriorating  stockpiles  from  two 
different  manufacturers,  and  so  on. 

We  distinguish  three  maintenance  models  in  which  only  minimal  repairs  are 
performed.  In  the  One-shot  Random  Repair  (OSR)  model,  a  coin  with  probability  X 
for  heads  is  tossed.  If  it  turns  up  heads,  a  unit  with  life  distribution  Gj  is 
placed  into  service  at  time  t  =  0  and  only  Gj-minimal  repairs  are  performed 
thereafter.  If  the  coins  turns  up  tails,  the  distribution  G^  is  used  in  place 


ely  Random  Repair  (CRR)  model,  a  coin  is  tossed  as  before 


of  Gj  in  the  above. 


and  depending  on  whether  it  turns  up  heads  or  tails,  the  unit  put  into  service 
at  time  t-0  has  life  distribution  or  G^.  At  each  new  failure,  a  coin  with 
probability  X  for  heads  is  tossed  independent  of  previous  history.  Depending 
on  whether  it  turns  up  heads  or  tails,  a  G^*miniinal  repair  or  a  G2"minimal  repair 
is  performed. 

In  the  Alternating  Repair  (AR)  model,  a  coin  is  tossed  as  in  the  previous 
models.  If  it  turns  up  heads,  then  a  unit  with  life  distribution  Gj  is  placed 
into  service  at  time  t=0.  From  then  on,  we  perform  minimal  repairs  in  alter¬ 
nate  order,  i.e.  G2-rainimal  repair,  Gj-minimal  repair,  G2-minimal  repair,  etc. 

If  the  coin  turns  up  tails,  the  roles  of  G^  and  G2  are  reversed. 

v;e  first  obtain  expressions  for  the  probability  generating  functions  of  the 
numbers  of  failures  before  t  in  the  above  three  models. 

Theorem  3.1.  Let  Nj(t)  be  the  number  of  failures  before  t  in  the  OSR  model. 

Its  probability  generating  function  is  given  by 

<{ij(t,  q)  =  X^(t)  +  (t),  where  p  =  1  -  q, 

and  its  expected  value  is 

(3.1)  ENj(t)  =  -  [X2,nGj(t)  +  UnG^2C^)]- 

Proof. 

This  follows  immediately  from  (2.3)  and  (2.6).  D 

For  the  CRP,  model,  we  first  show  that  the  probability  generating  function 
must  satisfy  certain  differential  equations.  We  assume  that  Gj  and  G2  are 


vv 
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absolutely  continuous  with  respective  densities  and  g^. 


Theorem  3.2.  Let  N^Ct)  be  the  number  of  failures  before  t  in  the  CRR  model. 
Its  probability  generating  function  is  given  by 

q)  =  XGj(t)A(t,  q) 

where  A(t,  q),  B(t,  q)  satisfy 

(3.2)  Gj(t)^  A(t.  q)  =  G^Ct)^  B(t.  q)  =  q[XA(t,  q)gj(t)  ♦  XBCt,  q)g2(t)] 

for  all  t  >  0, 


with  initial  conditions 


A(0.  q)  «B(0.  q)  »1, 


and  where  A  *  1  -  A . 


Proof. 

As  in  the  proof  of  Theorem  2.1,  we  let  be  the  time  of  the  n*^  failure 

and  let  F  be  the  distribution  of  X  .  Then 
n  n 


(3.3) 


P[N(t)=n]=F^^j(t)-F^(t) 


P(X  <tsx  ^,) 
^  n  n+l-^ 


(3.4) 


“  /JPCVl^tlX^.uXlF^tu) 


^  ^  /o 

°G  (u)  "  ®G  (u)  " 


Xa^(t)Gj(t)  ♦Xb^(t)G2(t) 


I.*,  /.  .v''  ‘ , '  V  %' 


where  a  (t)  ®  - , 

"  “GjCu) 

dF  (u) 

b  (t)  *  /  — - -  for  ns  1, 

"  °  G^Cu) 


and  a^Ct)  =bQ(t)  =  1. 

Differentiating  (3.3)  and  (3.4)  with  respect  to  t,  we  obtain  the  identity 


-dF^^l(t)*dFn(t) 

_  dF  (t)  _  dF  (t) 

=  XG  (t)^-2 - Xa„(t)g,(t)  4XG.(t)^2 - Xb^(t)g-(t) 

^  G  (t)  "  ^  2  G2(t)  "  2 


=  dF^(t)  -  Xa^(t)gj(t)  -  Xb^(t)g2(t). 


This  implies  that  for  nsi. 


(3.5) 


n*l 


=  /J 


Gj(u) 


(u)  dG  (u)  _  b  (u)  dG  (u) 

. . 

°  G  (u)  "  G  (u) 


and 


,a^(»)dGj(u)  b^MdCjCu) 
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G2(u) 


0  ;r 


G2(u) 


Let  A(t,  q)  =  I  q"a  (t)  and  B(t,q)  =  I  q”b  (t) . 
n=0  n=0 


Then  <|>(t,  q)  =  J  P[N(t)=n]q” 

n  =  0 


=  XG,(t)A(t,  q)  ♦XG,(t)B(t.  q)  . 


Mote  that 


a^(0)  *0  for  nil,  so  that 

A(0,  q)  »  B(0,  q)  = 1, 


A(t,  q)  =  1  +  q  I  q"a  (t) 
n  »  0 


»  dF  -  (u) 

r  _n  ft  n+1 


=  1M  I  q”/;-? 

n  =  0  G, 


Gj(u) 


by  (3.5). 


nP  r  ®T,f“5dG  (u) 

•i*q  I  q"  >/; 

n  =  0  "  Gj(u) 


b^(u)dG^(u)- 

®  Gj(u) 


Thus  Gj(t)“  A(t,  q)  =q[XA(t,  q)gj(t)  ♦X  B(t,  q)g2(t)l.  Similarly  we  have 

GjCt)^  B(t,  q)*q(XA(t,  q)gj(t)  ♦  X  B(t,  q)g2(t)].  Cl 

We  are  not  able  to  obtain  solutions  of  the  linear  homogeneous  partial  dif 
ferential  equations  given  in  (3.2)  for  arbitrary  G^  and  G2.  If  G^  and  G2  have 
proportional  hazard  rates,  i.e.. 


(3.6)  Gj(t)=G(t)  andG2(t)=G  (t)  for  some  0  <  6  <  1 , 

then  it  easy  to  solve  for  A(t,  q)  and  B(t,  q) .  The  solutions  are 


A(t,  q)  =6G”l(t)  ♦  8G“2(t) 


= 


where  a  ,  a_  are  the  roots  of  the  equation 


(a  qX)  (o  -  0  ♦  1)  ♦  aqXe  =  0 


and  6  =  1-$  is  chosen  such  that 

aj$  0^6 

Note  that  and  are  real  because 

[qX  +  (1  -  0)  ♦  qX0]^  -  4qX(l  -  0) 

=  [qX  -  (1  -  0)]^+  (qXe)^  ♦  2q)re[qa  ♦  (1  -  0)]  >0. 

Thus  the  probability  generating  function  of  N2(t)  is 

_  a.  $  . 

q)  =  txB+  ft) 

MXB* 

We  find  that  the  expected  number  of  failures  before  t,  under  the  assumption  (3.6), 
is  given  by 

—  2  — 

(3.7)  EN  (t)  =  -  3-^  tnG (t)  f  ^  [  1  -  G  (t) ] . 

^  X*XB  (X+X6)^ 

V/e  just  noted  that  we  were  unable  to  obtain  the  probability  generating  func¬ 
tion  of  N2(t)  for  arbitrary  and  G^.  However  we  are  able  to  obtain  EN2(t)  for 

arbitrary  Gj^  and  G^  in  Theorem  3.3  belo\v.  This  is  possible  because  the  identity 
(3.12)  (see  below)  allows  us  to  reduce  the  differential  equations  for  Aj(t)  and 
Bj(t)  which  determine  EN2(t)  (see  (3.15)  below)  and  solve  for  them.  Let  rj»*'2 
denote  the  respective  failure  rates  functions  of  G,  and  G_. 


Theoren  3.3.  In  the  CRR  nodel,  we  have 


(3.8)  EN2(t)  = 


-IXilnGj(t)  ♦  XtnG2(t)] 


fjGj(u)G2(u)[r  (u)  -r.,(u)]  .  Jp  tj;  ~x  - " 

°  ^  ^  ^  (iJ(v)G^(v) 


Proof. 


By  Theoren  3.2,  we  have 


(3.9) 


where 


q)  *  XGj(t)A(t.  q)  ♦  XG2(t)B(t,  q) 


(3.10)  A(t.  q)-l  =  qX/ 


^A(u,  q)dGj(u)  _  ^B(u,  q)dG2(u) 


Gj(u) 


Gj(u) 


B(t,  q)  -  1  *qX  /, 


^A(u,  q)dGj(u) 


G2(u) 


_  ^  B(u,  q)dG2(u) 

q^  Jp - z - 


G2(u) 


Let  A^(t)=A(t,  1)  and  Bj(t)=B(t,  1).  Then  (3.10)  implies  that 


Gj{t)A'(t)  =G2(t)BJ(t) 


(3.11) 


=  XAj(t),7j(t)  ♦  XEj(t)p.2(t) 


Putting  q =  1  in  (3.9),  we  have 


(3.12) 


XGj(t)Aj(t)  *XC2(t)Bj(t)  =  1. 


Solving  for  A, (t)  and  B, (t)  from  (3.9)  and  (3.10),  we  obtain 


Aj(t)Gj(t)C2(t)  =  XAj(t)[G2(t)Bj(t)  -  Gj(t)g2(t)]  +  g2(^)* 


which  implies  that 


G^(t) 


^  [Aj(t)Gj(t)G2\t)]  -  tA'(t)  ♦  X[r2(t)  -  rjCt)]Aj(t)} 


G^Ct) 


GjCt)  g2(t) 


GjW  Gj(t)G2(t) 


This  yields  the  solution 


Ai(t)=^ 


G2(t) 


G“(t) 


°  g[(u)G^*\u) 


Similarly 


^(t) 
B  (t)*:|: — 
^  ^(t) 


^  dGj(u) 


To  solve  for  EN-(t),  we  let  E(t)  =~  A(t,  q)  ]  ,  and  F(t)  =~  B(t,  q 

^  9c^  C|  s  1 

differentiating  (3.10)  with  respect  to  q  and  putting  q  =  l,  we  obtain 


A  (u)dG,(u)  _  B  (u)dG,(u) 


Gj(u) 


Gj(u) 


E(u)dG,(u)  _^F(u)dG2(u) 


S'-/  _  4. 


Gj(u) 


Gj(u) 


E(u)dG  (u)  _  F(u)dG,(u) 

=Ai(t)-i.x/;  2 


Gj(u) 


Gj(u) 


The  last  equality  follows  from  (3.10).  Thus, 


(3.13) 


E'(t)Gj(t)  =AJ(t)Gj(t)  ♦XE(t)gj(t)  *XF(t)g2(t). 


Similarly, 


(3.14)  F'(t)G2(t)  =  B'(t)G2(t)  ♦  XE(t)gj(t)  +  XF(t)g2(t) . 

Since  ^[EN2(t)]*^[^^2(^*  ‘»)lq=l] 


[XGi(t)E(t)  ♦XG2(t)F(t)l  by  (3.9) 


XE'(t)Gj(t)  ♦XF'(t)G2(t) 
-XE(t)gj(t)  -XF(t)g2(t) 


XA'(t)Gj^(t)  ♦XBj(t)G2(t) 


by  (3.13)  and  (3.14) 


=  XAj(t)gj(t)  +XBj(t)g2(t) 


by  (3.11) 


we  have 


(3.15) 


EN2(t)  =  X  /J  Aj(u}dGj(u)  ♦  X  /J  Bj(u)dG2(u) . 


Substituting  the  solution  of  Aj^(t)  and  3j(t),  we  obtain 


EN  (t)  =  X  / J  3:1 - dG  (u)  ♦  X  JJ  dG  (u) 

2  °gJ(u)  ^  ®^{u) 


.  ft  f“  "2^“’ 


-A  -T  dG  (v)dG  (u) 


Gj(u)  G^(v)G*^'‘(v) 

1 


dG,(v)dG  (u). 

rl^Xr.ftrXf  '  *  * 


^(u)  g}^^(v)G^(v) 


I 


Integrating  by  parts,  we  have 


_  t 

EN-(t)  *-tnG,(t)  -  X/;  -  dG  (u) 

2  1  0  G  (u)  ^ 


(3.16) 


-,t  "iM 

♦X/^  -  dG,(u) 

°  G  (u)  2 


Gj(u) 


♦  ^  /J  Z= -  82^“^  ‘^X - U  -\'*T  -X - 

°  G^(u)  2  gJ(u)  ^  °  g}*^(v)G^(v) 


dG,  (V) 


Interchanging  Gj  with  G2  and  X  with  X  in  (3.16),  we  obtain  the  following  alter¬ 
native  expression  for  EN2(t): 

EN-(t)  =  -)lnG-(t)  -  X  /J  ^ -  dG,(u) 


'»  G,(u, 


(3.17) 


rt 

*xr  ^ -  dG  (u) 

°  G  (u)  ^ 


Gi(u) 


G^u) 


ft  ''tv-/  u  — 2 


g;(u) 


dG^(v) 
u  2 


G^(u)  ^  ”  G^(v)Gp(v) 


The  expression  for  EN2(t)  in  (3.8)  is  obtained  from  (3.16)  and  (3.17)  by  appro¬ 
priate  averaging.  0 

We  next  present  analogous  results  for  the  AR  model.  Let  Nj(t)  be  the  num¬ 


ber  of  failures  before  t  in  the  AR  model. 


mmmmm 


Lemma  3.4.  Consider  the  conditional  distribution  of  given  that  a  head 

turns  up  on  the  first  toss,  which  is  the  same  as  saying  that  the  first  unit  has 
life  distribution  Its  probability  generating  function  is 

Gj(t)A(t.  q)  ♦G2(t)B(t.  q)  , 
where  A(t,  q)  and  B(t,  q)  satisfy 

q)=q8(t,  q)g2(t), 

B(t,  q)=qA(t,  q)gj(t), 

with  initial  conditions 

A(0,  q)  =  1,  B(0,  q)  =0. 

The  expected  value  of  N,(t)  given  that  the  first  unit  has  life  distribution  G,  is 

1 

\ 

-tnG^Ct)  ♦/J  Gj(u)G2(u)[rj(u)  -r2(u)].  l  +  /“  - du. 

G,(v)G2(v) 

Proof. 

The  proof  is  omitted  because  it  is  similar  to  the  proof  of  Theorems  3.2 
and  3.3.  □ 

The  expected  value  of  N^Ct)  follows  immediately  from  Lemma  3.4.  We  state 
this  result  without  giving  the  proof. 
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Theorem  3.5.  In  the  AR  model. 


(3.18) 


ENj(t)  =  -[XtnG^Ct)  ♦  XtnGjCt)] 


rt  -  -  f  _  dG  (V) 

-/q  Gj(u)G2(u)[rj(u)  -r2(u)]]-(X-X) 

Gj(v)G2( 


-  ^ /qZ - =5 - M'"* 

°Gj(v)g2(v) 


We  may  now  compare  the  expected  number  of  failures  among  the  three  models 
for  the  case  where  A  =  X  =  l/2.  We  further  assume  that  the  failure  rate  functions 
satisfy 


rj(t)2:r2(t)  for  all  t^O. 

Then  it  follows  from  (3.1),  (3.8)  and  (3.18)  that  the  expected  number  of  failures 
before  t  is  the  largest  in  the  OSR  model.  To  compare  the  AR  model  and  the  CRR 
model,  it  is  clear  from  (3.8)  and  (3.18)  that  the  expected  number  of  failures 
before  t  is  smaller  in  the  AR  model  than  in  the  CRR  model  if  and  only  if 


(3.19) 


,  rt?-l/2,  ,-1/2,  ,,  ,,  ,  ,,  ru  1^  ^  ^2^  ^ 

1/2/  G  (u)G  '  (u)(r  (u)  -  r  (u)]- 

i  1  ^  [  '^Gj'^(v)G2'  (\ 

s  j:  G,(u)G  (u)[r  (u) -r_(u)]{j«-i - ^ - dv  du. 

^  ^  Z  0  Gj(v)G2(v) 


dv  >  du 


We  were  not  able  to  prove  (3.19)  for  all  tiO.  However,  for  small  values  of  t, 
namely  when  Gj(t)G2(t)  s 1/4,  it  is  straight  forward  to  show  that  (3.19)  is  true 
for  arbitrary  G^  and  G2.  Therefore,  if  the  process  is  to  be  run  for  a  short 
time,  AR  model  is  preferable  to  the  CRR  model.  Under  the  proportional  hazard 


asstimption,  namely  (3.6),  the  expected  number  of  failures  before  t  is  smaller 
in  the  AR  model  than  in  the  CRR  model  for  all  t,  since  (3.19)  can  be  reduced  to 
the  inequality 

1/2(1  -  G^^®(t)]  i  1  -  G^^'^®^^^(t), 
which  holds  for  all  t^O. 

The  BP  model  cf  Section  2  can  be  extended  to  include  perfect  repair  and 
two  kinds  of  minimal  repairs.  Once  again  epoch  of  perfect  repairs  form  regen¬ 
eration  points.  The  tail  of  the  distribution  of  the  waiting  time  between  per¬ 
fect  repairs  is  again  the  probability  generating  function  of  the  number  of 
failures  in  the  same  repair  model  when  perfect  repairs  are  excluded. 

Repair  models  with  k  minimal  repairs  where  k  s  3  can  be  obtained  by  straight 
forward  extensions.  Differential  equations  for  the  various  quantities  involved 
similar  to  those  found  in  Theorem  3.2  are  easily  obtained.  It  turns  out  that 
these  linear  homogeneous  differential  equations  are  not  easy  to  solve  even  in 
the  proportional  hazard  case. 
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